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Abstract 



As a response to ligand binding, transmembrane cell receptors often en- 
hance their clustering, or oligomerization, during the signalling process. Here 
we present a statistical mechanical model which combines the aspects of clus- 
tering and signalling. In this model, receptors float on the surface, while 
for two neighboring receptors, there is an interaction energy dependent on 
their conformational states. On the other hand, ligand binding of a receptor 
shifts the energy difference between the two conformational states. Due to 
thermal fluctuation, the effects of clustering and signalling are statistical av- 
erage quantities. This model reduces to a floating Ising model with a random 
field. We calculate the signalling in a grand canonical ensemble mean field 
approach, using Hubbard-Stratonovich transformation and replica method. 
Monte Carlo simulations are also performed. Essential biological features are 
obtained in our model. 
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1. INTRODUCTION 



On the surface of a cell, there are many transmembrane receptor molecules. They are 
"sensory machines" of the cell. Binding extracellular ligands changes the conformational 
states or, in our point of view, the probability distribution of the conformational states. 
Consequently a cascade of responses occurs inside the cell. This is called signalling, which 
is often quite sensitive. An example is bacterial chemotaxis Q: through the detection by 
its receptors, a bacterium swims towards higher concentration of a chemical it likes or a 
lower concentration of a chemical it dislikes. In this case, thousands of receptors cluster 
together at a pole of the cell. Recently, based on detailed analyses on the experimental 
findings, I proposed that cooperative receptor-receptor interaction enhances the signalling, 
and constructed a statistical mechanical theory for this process [0-0]. My idea has since 
received considerable attention []. In [|]], good consistency was found in comparing this 
theory with experimental results [|5|-|TU|, and it was also suggested what kinds of experiment 
are needed. 

In the case of many other cell receptors, however, usually they float on the cell membrane. 
A ubiquitous phenomenon is that ligand binding enhances the clustering of receptors during 
the signalling process. Investigations on this issue are widely reported for various receptors, 
for example, integins, which attach cells to the extracellular matrix, tyrosine kinase receptors, 
cytoskine receptors, growth factor receptors, tumor necrosis factor receptors, including death 
receptors, antigen receptors such as T-cell receptors, bone morphogenetic protein receptor, 



G-protein-coupled receptors, etc pT|-[19|. In our previous works @-§J], the position of each 
receptor is fixed, and there is an a priori cluster. In the present article, we shall consider 
statistical mechanics of floating receptors. Based on our model, it is proposed that signalling 
sensitivity and clustering enhancement, in certain cases, are likely two consequences of the 
same mechanism based on conformation-dependent local receptor-receptor interaction. 

In Sec. 2, we give a detailed description of the model. A mean field calculation of the 
signalling, as a statistical average quantity quantity, is given in Sec. 3. In Sec. 4, we discuss 
receptor clustering. Results of Monte Carlo simulations are given in Sec. 5. Finally a 
summary and discussion is made in Sec. 6. 



2. THE MODEL 

Consider a collection of receptors on a lattice. Fig. 1 is an illustrative configuration. 
Let the coordinate number be is, which is 6 for a honeycomb lattice and is 4 for a square 
lattice. Suppose that the receptors can move around on the lattice. Hence at an instant, 
a receptor has < v nearest-neighboring receptors. The conformational state of the receptor 
at lattice point i is represented as Vi, the value of which is either V° or V 1 , if we follow 
the two-state convention. ]/] can be understood as the position of the receptor molecule. 



The reference was mostly made to a publication which popularizes some rudiments of my idea 
[T. Duke and D. Bray, Proc. Nat. Acad. Sci. 96, 10104 (1999)], and the attribution of credit was 
inappropriately made in T. Duke, N. Le Novere and D. Bray, J. Mol. Bio. 308, 541 (2001). 
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Obviously, V, is influenced by whether this receptor is bound to ligand. We suggest that 
Vi is also influenced by the conformational states of its nearest-neighboring receptors, due 
to receptor-receptor interaction which is dependent on their conformational states. We use 
Hi to represent whether there is a ligand binding, that is, H { = H ^ if the receptor is 
bound to ligand, otherwise Hi = 0. H can be understood as the force or torque generated 
by the ligand binding. HVi is the energy due to ligand binding. Given the ambient chemical 
concentration [L], the occupancy, i.e. the fraction of receptors with ligand bound, is 

- [L] (D 



where is the dissociation constant. Therefore each Hi is randomly distributed between 
H and 0, with probabilities c and 1 — c, respectively. Formally, the probability distribution 
can be written as 

p{Hi) = c5(Hi -H) + (l- c)5(Hi), (2) 

where 5(x) = if x ^ 0, while J 5(x)dx = 1. Since usually K& < [L], it is reasonable to 
regard the randomness of Hi as quenched 0,|J. 

We associate each site i with a variable Hi, which is 1 if there is a receptor on this site, 
and is otherwise. In reality, of course Vi and Hi are only defined when there is a receptor 
on site i. However, because the presence or not of a receptor at site i has been determined 
by Hi, for the convenience of treatment, equivalently we may associate Vi and Hi to each site 
i, no matter whether there is a receptor on site i. Since the receptors are allowed to diffuse, 
rij is time-dependent. One may constrain that the total number of receptors, N = J2i n i, is 
conserved. It is also possible that the receptors may get out of the relevant region, hence N is 
not conserved. In any case, we may use grand canonical ensemble approach characterized by 
the chemical potential, which, in the first case, is given by the constraint on the conservation 
of receptor number. 

Therefore, in absence of thermal fluctuation, we have 

Vi = + Hi- Ui), with^(x) = { j| ^ < o > (3) 

j < — 

where the summation is over the nearest neighbors j of site i, Ui is a threshold value, is a 
coefficient for the coupling between nearest neighbors. J2j tij n jVj is the influence, i.e. force 
or torque, on the receptor at site i from the nearest-neighboring sites. Hi is the influence, 
i.e. force or torque, generated by the ligand binding of the receptor at % itself. Eq. (^) means 
that if the total influence on a receptor, as the sum of that due to its ligand binding and 
that due to the interaction with the nearest-neighbors, is larger than a certain value Ui, the 
receptors is in state V 1 . Otherwise it is V°. Equivalently, with ty = tji and ta = 0, the 
dynamics is determined in the following Hamiltonian, 

H = runjViVj - H i n iV + U ]T mVi, (4) 

(ij) i i 

where (ij) represents pairs of nearest neighbors, £y and U have been set to be constants. 
Eq. (fD gives the total energy of the system, with the kinetic energy neglected. The first 
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term is the total interaction energy of nearest-neighboring pairs. The second term is the 
energy due to ligand binding. The third term is the original energy, which, together with 
the first term, determines the probability distribution of the conformational states of all 
the receptors in the absence of ligand binding of any receptor. Without ligand binding and 
coupling with others, the energy of a receptor is UV° for state V°, and is UV 1 for state V . 
With ligand binding but without coupling, it is (U — H)V° for state V°, and is (U — H^ 1 
for state V 1 . The interaction energy of the pair (ij), if sites i and j are both occupied by 
receptors, is —tViVj, the value of which is dependent on both Vi and Vj. Since ligand binding 
can cause more Vi to flip from V° to V 1 , it is required that (H - U)^ 1 - V°) > 0. Precise 
form of TC is subject to future experiments. However, (H) captures the essential elements 
while allows a theoretical treatment. 

For convenience, we transform to a spin language, with the spin defined as 

Hence it is Si = 1 if Vi = V 1 , and is Si = — 1 if V, = V°. In order to simplify the mathematics 
without loss of the essence, suppose V 1 = —V°. Therefore we obtain the following simplified 
Hamiltonian 

H = -J niUjSiSj - Y BiUiSi, (6) 

<ij> i 

with J = t(V rl -1/°) 2 /4, B t = (H l -U)(V 1 -V°)/2, which is thus B x = (H - U^V 1 - V°)/2 
with probability c and is B 2 = — UiV 1 — V°)/2 with probability 1 — c, i.e. 

p(Bi) = c5(B t - B x ) + (1 - 0)5(3, - B 2 ). (7) 

We call this model, defined by (0) and (0), floating Ising model with a quenched random field 
at each site. Note that it is not a diluted Ising model, since rii is not a quenched disorder, 
but is a dynamical variable on the same footing as Si. 

It should be mentioned that a cell signalling process often involves different and largely 
separated time scales. On a much longer time scale the degrees of signalling and clustering 
evolve after the stimulus, till they return to the pre-stimulus levels. So there is a piece of 
short-time scale physics, and there is a piece of long-time scale physics. Here it suffices to 
discuss the the short-time scale physics, i.e. the quasi-equilibrium part of the theory. With 
an incorporation of the feedback, the long-time scale behavior is similar to that described in 
Ref. . With the adaptation process taken into account, the ligand concentration above 
should be replaced as the change of the ligand concentration. However, as studied in H|J, 
perfect adaptation is equivalent to no ligand binding. 

Guo and Levine studied the clustering of receptors under the assumption that the in- 
teraction between neighbors depends on whether the receptors are liganded [P0[] , while the 
signalling problem was not considered. As described in our model, We think that it is 
more likely that the interaction depends on the conformations of the receptors, rather than 
whether the receptors are liganded. In addition to the consistency with the situation of 
chemotaxis JMJ], this is also supported by the strong experimental evidence that oligomer- 
ization is stabilized by receptor-receptor interaction promoted by conformational changes of 
the receptors 0-|l|,p|,|I| . 
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3. MEAN-FIELD CALCULATION OF SIGNALLING 



The activity of the system is (Yli n iVi)i where the summation is over all the lattice 
sites, (. . .) represents ensemble average, i.e. average over the probability distribution of all 
possible configurations of the system. It is straightforward to write it as (V 1 — V°)W/2, 
with W = (J2i n>iSi) = Im, where m = (riiSi) is the average spin per site, I is the number of 
lattice sites. The signalling can be characterized by the change of W when the occupancy 
is changed from c = to c, i.e. W(c) — W(c = 0) = I[m(c) — m(c = 0)]. 

Therefore we set out to calculate m as a function of c, using grand canonical ensemble 
approach. The grand partition function is 

* = IE E E xn^^e-^, (8) 

i m=0,l 5j=-l,l i 

where \i is the chemical potential, Ti is as given in (^j). Using Hubbard-Stratonovich trans- 
formation, we obtain 



A /(II d</>i) exp[-i £ frK^h + + PBJniSi] 

i i,j i 

where A = {2Tr)- I / 2 {detK)' 1 ' 2 . K = /3(Jy) is a matrix, 1 is its inverse. Consequently 

= A fill dfa) exp(-~ £ ^ Vi) H 2 + ePi**-** + e^"^]. (9) 

i i,j i 

Based on In J? = Hm n ^o(Z n — the replica method yields 

InZm}) = lim - U n J ([] #?) exp(-^ 2 IIP + e^+^ + eP^tf-PBi] - 1 



(10) 



where a = 1, • • • , n denote the replicas. 

Averaging In Z over the Bi distribution (|7|), one obtains 

hTZ= [\nZ({B l })Hp(B l )l[dB l 

J i i 

= lim -{A n f (II^?)«p(-i £ Wffi 

n— »u 77, j ^ 



n 



Cl rj(2 + e /3M+(^+/3Bi) + ^n-wz+m)} _|_ C2 rj(2 + eft*+-w?+/»£a) + e /?/i-w+/»sah 



where, for notation convenience, we have set Ci = c, and C2 = 1 — c. 

The mean field value of 4>f can be obtained by deciding the saddle point of the integrand 
in Eq. (|TT|) , for which one can obtain 



v-^ K-ifo c ifi 1 9i + c 2 / 2 n x gi no , 

SA -*" cJt + cfi ' <12) 
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where $ n is the saddle point value of (ftf , corresponding to n replicas, fk = 2 + e /3M+ (* n+/3Bfc ) + 
e /9^-(* n +/9B fc ) > g k = e ^+{*n+PB k ) _ e ^-(^ n +pB k )^ ( k = Taking n -f limit, we have 

Ej Kjj 1 ® = c 1 g 1 /J l + c 2 ^ 2 /7 2 , where $ = lim n ^ $n, 7 fc = 2 + e /^+(*+/^) + ^-^ B "\ 
g k = e/?Ai+(*+/3B») _ e ^-(*+^B*) j (fc = i >2 ). Therefore 

$ = ^J^ + c 2 ^. (13) 

is the mean-field value of the Hubbard-Stratonovich field. The saddle point approximation 
implies that 

= lim ^ |exp(-| £ $ ^ )(ci7? + cjj)' - 1 J . (14) 

m, as a measure of the activity, is calculated in the following way. For given {Bi}, 
the ensemble average of Si is (Si){ Bi } = dlnZ({Bi})/d(/3Bi) = dlnZ '({Bi}) / d(<f>i) = 
J2j Kij 1 (4>i) {B t } , whose average over p({{B { }) gives the m. Thus m = (Si) {Bz} = Y^jKOj 1 ®, 
where $ is the mean-field value of the Hubbard-Stratonovich field, as given in ([131) . Therefore 
$ = f3vJm. Eq. fllTf ) can thus be written as 

e /3[t+f}vJm+pB k _ e fj^L-fjuJm-l3B k 
m ~ ° fc 2 + e Pn+PvJm+pB k _|_ e P^-PuJm-l3B k ' (^) 



The average receptor number is (N) = <91n Z/d((3/j,), and can be evaluated as 

gf3fi+f3vJm+f3B k _|_ ^/3fi—/3uJm—/3B k 

(16) 

( |15D and ([IB]) can also be obtained in a simple mean field approach regarding H — fiN 
as an effective Hamiltonian, and considering four possibilities of the state (n^, Si) at site i, 
subject to an effective field uJm + Bi. 

Hence we have obtained an analytical expression for the signalling, as a statistical average 
quantity. If the receptor number is conserved in the system, (|T6|) gives /iasa function of 
the receptor number. 

Sensitivity of the signalling to ligand concentration is given by dm/dc, which can be 
arbitrarily large. 

For annealed randomness, which may apply to the rare case in which [L] is comparable 
to Kd, one may obtain 

Y,k-1 2 Ck(e l3 ^ +l3uJm+l3Bk — e P^-PvJm-f3B k ^ 
m ~ J2k=l 2 Cfc(2 + e P»+PvJ™+/3B k _|_ e /3fi-l3uJ m -l3B k ^ ' ( 17 ) 



and 



V r, tpPv+PvJm+pBy. _|_ P/j,-j3uJm-j3B k \ 

/ N \ _ j 2^=1,2 Cfcje +e )_ 

In the high temperature limit, the difference between flTj ) and (p|), and that between 



fllq) and (|T8[) tend to diminish. 
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4. CLUSTERING 



The conformation- dependent interaction not only enhances signalling sensitivity, it is 
also responsible for clustering. The first term of the Hamiltonian implies that, in order 
to minimize the Hamiltonian, the receptors tend to aggregate together to maximize the 
number of nonzero riiUj with Si = Sj for neighboring (ij). Therefore receptors with the 
same conformational state tend to cluster. However, because of thermal fluctuation, they 
cannot all cluster together, since there is a probability distribution for various configurations. 
As an illustration of the situation, a snapshot obtained in a Monte Carlo simulation is shown 
in Fig. 1. 

The second term in (^) determines how ligand binding affects the clustering situation. It 
can be seen that U determines the bias of the distribution of the receptor state {Si} between 
1 and —1 when there is no ligand binding (Hi = 0). That is, if U = 0, for not-too-strong 
interaction compared with the temperature (/3J)Q, there is an equal number of receptors 
with spin 1 and those with —1. If U > 0, there are more receptors with spin — 1. U also 
determines whether ligand binding enhances or suppresses the clustering. For example, we 
may consider two typical cases: 

Case (i),U = H . According to fl3|), this means that merely a ligand binding of a receptor 
is enough to change its conformational state. In this case, B\ — 0, B 2 — —HiV 1 — V°)/2. 

Case (ii), U = 0. According to (|4]), this means that in the absence of ligand binding 
(Hi = 0), there is no bias in the state distribution. In this case, B\ = H(V l — V°)/2, B 2 = 0. 

These two cases can map to each other with the transformation H — > —H and c — > 1 — c. 

In case (i), since B 2 < 0, in the absence of ligand, all Si tend to be close to —1. Ligand 
binding causes more Si to be 1. Because receptors with the same value of Si tend to 
cluster, those with different Si tend to be disconnected, the conclusion is that ligand binding 
suppresses clustering. 

In case (ii), if there is no ligand binding, there is no biasing field, so there are equal 
probabilities for Si to be 1 and —1, consequently the clustering is minimized. When there 
is ligand binding, the clustering is enhanced. 

For other values of U, i.e. U ^ and U ^ H, we have the following general picture. 
Without ligand binding, there is a uniform field B 2 at every site. In the presence of ligand 
binding, there is a random field as given by (0). With a rough mean field estimation, one 
may see that whether the clustering is enhanced depends on whether (1 — c)B 2 + cB\ has 
larger absolute value than \B 2 \. Thus only if H > 2U /c, is the clustering enhanced by ligand 
binding. Case (ii) belongs to this region. If U — H/2, we have B\ = —B 2 . 

The clustering of receptors can be studied quantitatively by defining a clustering corre- 
lation function as 

C(r) = (^+;>-(n> 2 
(n) 2 

where the average is over different sites % and different directions of r first, and then over 
the thermodynamic ensemble, i.e. different possible configurations at the same temperature. 



2 There is no such a limitation if the adaptation on a long time scale is considered, see Ref. HQ. 
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(n) = (N)/I is the density of receptors on the lattice. If there is no clustering correlation, 
(riini +r ) = (n) 2 , consequently C(r) = 0. This correlation function measures the deviation 
from the non-correlating case and allows comparison of situations for different receptor 
densities. This definition is similar to the two-point correlation function in the study of 
galaxy clustering. One may also define higher-order correlation functions. We leave for 
future investigations the analytical calculation of the clustering correlation function, and 
turn to numerical simulations here. 

5. RESULTS OF MONTE CARLO SIMULATIONS 

To investigate the extent of clustering and to calculate the activity, we have done Monte 
Carlo simulations using the Metropolis algorithm on a square lattice. In the simulations, 
we conserve the receptor numbers. We have specifically studied Case (ii), from which one 
may obtain the results for case (i) simply by changing c to 1 — c. The results show strong 
correlations for small values of r and weak correlations for large values of r. The decay looks 
like exponential, as expected for a non-critical system. 

We studied correlation functions for different values of the ligand binding fraction c, with 
a same receptor density (n). See Figs. 2(a) and 2(b) for results under different values of the 
coupling- noise-ratio f3J. From the correlation function for small values of r, it is clear that 
the larger the ligand binding fraction c, the larger the correlation. This confirms the above 
analyses. 

Comparing Figs. 2(a) and 2(b), it can be seen that with larger coupling-noise ratio (3 J, 
as in Fig 2(a), the correlation function C(r) is larger at small values of r, while smaller at 
large values of r, indicating that clustering is stronger. On the other hand, for larger (3 J, 
the decay of C(r) with r is faster, because it is further from criticality. 

We also studied the correlation functions for different values of receptor densities (n) with 
a same ligand binding fraction c (see Fig. 2(c)). It is shown that the smaller the receptor 
density, the larger the correlation function. This can be understood, since the larger the 
density, the less freedom two receptors can approach each other. 

In the simulations, we also obtained the activity. The activity W can also be written as 
Y!i Si, where J2'i represents summation over the receptors instead of the lattice sites. Hence 
it equals NM, where M is the average activity per receptor. 

Fig. 3 gives the relation between the average activity per receptor M and the ligand 
binding occupancy c, for different values of coupling- noise ratio f3J . From the plots, we see 
that the activity increases with the ligand binding fraction c, with (3 J, and with the receptor 
density (n), as consistent with the mean field solution. 

6. SUMMARY AND DISCUSSION 

In this paper, we propose a statistical mechanical theory which accounts for clustering 
and signalling of receptors in a same framework. Clustering and signalling of a network of 
receptors are treated as statistical average quantities. In our model, which can be reduced 
to a so-called floating Ising model in a random field, the interaction energy between neigh- 
boring receptors depends on their conformational states, therefore through this interaction, 
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the conformational state of one receptor influences those of its neighbors. On the other 
hand, since the receptors are allowed to move on the membrane, the receptors with a same 
conformational state tend to cluster together, in order to decrease the total energy of the 
system. Therefore clustering and signalling are unified as two consequences of the same 
coupling between receptors. 

According to our theory, clustering exists even in absence of the ligand binding. This 
is consistent with a recent experimental finding [I8|. We have studied clustering, based on 
an appropriate definition of clustering correlation function. Monte Carlo simulations were 
made. To obtain the situation that ligand binding enhances the clustering, a parameter, 
namely the threshold value U for the change of the conformational state, which determines 
the state distribution in the absence of ligand binding, must be within a certain range. If the 
values of the variable characterizing the two conformational states are symmetric, i.e. with 
a same magnitude and with opposite signs, then a simple possibility is that this threshold 
value is zero. We note that it was found experimentally that, in the absence of ligands, the 
receptors are hindered to cluster by certain inhibitors, which are squeezed out when ligands 
bind |21[ . In such a case, it may be constrained that clustering is always enhanced by ligand 
binding. 

By using a mean field theory based on Hubbard- Stratonovich transformation and replica 
method, as well as by using Monte Carlo simulations, we also studied the activity as a 
function of ligand binding fraction c. 

This model is also be interesting from the point of the view of statistical mechanics. 
On the other hand, it is straightforward to make appropriate extensions of this model, for 
example, to put in more details of the realistic systems. The adaptation can be studied 
by straightforwardly generalizing our previous approach based on a counteracting field as 
a feedback from the signal to the field ]3[]. In this case, both signalling and clustering are 
adapted through a feedback, on a long time scale. 

A simple and direct experimental test of our theory is to examine whether conformational 
state of a receptor can be changed by ligand binding of its nearest-neighboring receptors. It is 
interesting to study the forces generated by ligand bindings of a receptor and its neighbors. 
Parameters in our model need to be measured. The techniques used in studying single 
molecules may be useful for this subject. 
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Figure Captions: 

Fig 1. A configuratiorial snapshot of 300 receptors on a 50 x 50 lattice. It is in 10000 
Monte Carlo steps after an initial random configuration. An empty circle represents a 
receptor with Si = —1 and Bi = 0, a filled circle represents a receptor with Si = —1 and 
Bi = B, an empty square represents a receptor with Si = 1 and Bi = 0, a filled square 
represents a receptor with Si = 1 and i?j = £>. The probability for Bi to be 5 is 0.8. 
j3B = 2, /5J = 0.8. Obviously, increasing (3 J enhances clustering. 

Fig 2. Monte Carlo calculation of the clustering correlation function C(r) for receptors 
on a 50 x 50 lattice. The result is obtained by averaging over different sites with the 
same distance, and over 1000 Monte Carlo steps following 10000 initial steps to approach 
equilibrium. (3B = 2. (a) f3J = 1.6, there are 300 receptors, results for different values 
of ligand binding fraction are displayed as different symbols, (b) f3J = 0.8, there are 300 
receptors, results for different values of ligand binding fraction are displayed as different 
symbols, (c) [3J = 0.8, the ligand binding fraction c = 0.8, results for different total 
numbers of receptors are displayed as different symbols. 

Fig3. Monte Carlo calculation of the activity of the receptors on a 50 x 50 lattice. 
For comparison between different densities of receptors, we give the relation between the 
average "spin" per receptor M and the ligand binding occupancy c. The result is obtained 
by averaging over all receptors and over 1000 Monte Carlo steps following 10000 initial 
steps to approach equilibrium. f3B = 2. (i) iV = 300, (3 J = (•), i.e. there is no receptor- 
receptor interaction, (ii) N = 300, [3 J = 0.4 (+). (iii) N = 300, f3J = 0.8 (x). (iv) N = 300, 
f3J = 1.6 (©). In this case, the deviation of M(c = 0) from is spontaneous "magnetization" 
due to strong coupling, but it is always brought to by the negative feedback on a long time 
scale |Il . (v) N = 900, f3J = 0.8 (<>). 
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